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ABSTRACT 


A  further  clarification  of  the  volume  from  which  data  originate  in 
the  dual  scatter  laser  Doppler  velocimeter  is  presented.  There  are 
several  factors  that  specify  the  primary  dimensions  of  the  probe  volume 
which  in  some  cases  are  totally  specified  by  the  collector  optics.  It 
was  found  experimentally  that  the  lens  aperture  of  the  light  collecting 
optics  has  a  primary  effect  on  the  volume  from  which  data  originate. 

An  equation  is  derived  for  reducing  this  data  volume  for  both  the  on- 
axis  and  off-axis  cases.  Experimental  verification  is  presented  and 
compared  to  the  derived  equation. 
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SECTION  I 
INTRODUCTION 


The  dual  scatter  laser  Doppler  velocimeter  (LDV)  has  been  routinely 
used  to  acquire  velocity  information  from  extremely  small  spatial  regions 
in  flow  fields.  Previous  attempts  to  accurately  obtain  a  correlation  be¬ 
tween  the  experimentally  deduced  volume  and  the  calculated  volume  have 
been  exceedingly  difficult.  The  experimental  values  obtained  for  probe 
volumes  were  found  to  be  much  smaller  than  the  predicted  values.  This 
is  particularly  true  in  the  case  of  the  dual  scatter,  backscatter,  LDV 
system.  In  this  particular  case,  a  concentric  light  collection  system  is 
used  in  conjunction  with  the  laser  beam  focusing  system  (on-axis  back¬ 
scatter  unit).  The  important  aspect  is  that  there  are  two  distinctively 
separate  spatial  volumes  under  consideration.  The  volume  defined  by  the 
two  laser  cross  beams  is  termed  the  geometric  volume;  the  volume  asso¬ 
ciated  with  the  light  collecting  optics  is  termed  data  volume.  Henceforth, 
when  reference  is  made  to  the  probe  volume,  it  is  to  mean  the  effective 
volume  from  which  information  is  extracted  with  a  given  optical  configur¬ 
ation.  In  the  following  discussions,  all  optical  ray  boundaries  are  defined 
as  the  1/e^  (0.  135335. . .)  amplitude  having  assumed  a  Gaussian  curve  as 
the  amplitude  distribution. 


SECTION  II 

THEORETICAL  DETERMINATION  OF  PROBE  VOLUME  DIMENSIONS 


In  general,  probe  volumes  have  been  specified  as  the  maximum 
radial  diameter  and  the  maximum  axial  length  of  the  geometric  probe 
volume  defined  below.  This  definition  can  be  misleading  where  probe 
volumes  are  to  be  minimized  when  only  the  reduction  of  the  geometric 
volume  is  considered.  It  is  important  not  to  overlook  the  reduction  of 
the  probe  volume  that  is  also  obtained  with  the  light  collecting  optics. 
This  feature  is  emphatically  demonstrated  in  Fig.  1  (Appendix  I)  which 
contains  characteristic  curves  of  the  axial  dimension  of  the  probe  volume 
obtained  from  an  on- axis  and  off-axis  optical  arrangement  coupled  with 
an  0 . 0 1 8  -in.  aperture . 


2.1  GEOMETRIC  PROBE  VOLUME 

The  geometric  volume  is  formed  by  the  intersection  of  the  two  laser 
cross  beams.  Each  laser  beam  is  focused  to  its  minimum  value  at  the 
region  of  intersection.  Hence,  since  each  laser  beam  may  be  described 
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by  a  Gaussian  intensity  distribution,  the  focused  region  may  also  be 
described  by  a  composite  Gaussian  intensity  distribution.  The  minimum 
focal  diameter  is  determined  for  a  diffraction-limited  focusing  system 
by 


2b  ■  £k  A 

0  rr  2b 


<1) 


where  b0  is  the  radius  of  the  focused  spot,  fc  the  focal  length  of  the 
focusing  lens,  X  the  laser  radiation  wavelength,  and  2b  the  laser  beam 
diameter  at  the  focusing  lens.  It  has  been  shown  (Ref.  1)  that  three 
dimensions  are  associated  with  the  laser  beam  crossover  region  shown 
in  Fig.  2.  These  dimensions  describe  an  ellipsoid  at  the  1/e2  inten¬ 
sity  points.  The  determination  of  the  Ax,  Ay,  and  Az  dimensions  may 
be  easily  deduced  (Fig.  2)  as 

Ax  =  2b0 

Ay  =  — ^ —  /2\ 

cos  9/2 


Az 


2b0 

sin  9/2 


where  9  is  the  angle  between  the  two  cross  beams.  Nomograms  have 
been  developed  (Ref.  4)  that  allow  rapid  determination  of  these  geometric 
probe  volume  parameters. 


2.2  DATA  PROBE  VOLUME 


All  of  the  velocity  information  is  contained  totally  within  the  geo¬ 
metric  volume;  however,  the  data  volume  is  only  a  portion  of  the  entire 
geometric  volume  from  which  the  velocity  information  is  extracted. 


In  attempting  to  minimize  the  geometric  volume,  the  major  problem 
lies  with  the  resolution  of  the  axial  dimension,  or  the  depth  of  field, 
which  is  much  larger  than  the  radial  dimensions;  specifically,  from 
Eq.  (2) 


Al  =  2b'°  /UJ\  =  CSC  9/2 

Ax  sin  9/2  \2b0/ 


(3) 


i 


Equation  (3)  indicates  how  much  greater  this  dimension  is  compared  to 
the  (effective)  diameter  of  the  ellipsoid.  A  velocimeter  with  a  3-deg 
angle  (typical)  between  beams  has  an  axial  dimension,  or  depth  of  field, 
38  times  greater  than  the  effective  diameter.  The  lens-aperture 
arrangement  of  the  collector  optics  is  one  technique  that  may  be  used  to 
reduce  this  axial  dimension. 
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2.3  PROBE  VOLUME  DETERMINATION 


From  past  experience,  background  light  has  a  considerable  effect 
on  the  signal -to -noise  ratio  of  the  detected  Doppler -shifted  frequencies. 
Even  though  efforts  are  generally  made  to  reduce  overhead  lighting, 
some  experiments  have  been  set  up  with  no  consideration  to  overhead 
lighting  effects.  It  was  found  that  with  an  adequate  light  shield,  and  a 
black  curtain  placed  in  a  suitable  location,  no  detectable  light  could  be 
measured.  Figure  3  is  the  schematic  of  this  setup.  This  experiment 
suggested  that  only  light  emanating  from  the  centerline  of  the  lens 
arrangement  will  contribute  to  the  noise  term  of  the  velocimeter. 
Therefore,  the  effect  of  unfocused  light  on  the  centerline  of  the  lens 
arrangement  was  studied  also. 


2.3.1  Theoretical  Analysis,  Axial  Dimension 


The  diameter  of  the  spot  of  a  focused  laser  beam  may  be  deter¬ 
mined  from  Maxwell  equations.  The  laser  beam  propagates,  say,  in 
the  +z  direction  and  spreads  only  by  diffraction.  The  local  direction  of 
propagation  (Ref.  2)  is  hyperbolic  in  nature,  and  over  a  spherical  sur¬ 
face,  of  constant  plane,  the  field  amplitude  is  given  by 

Em,n  =  E„  IIm(vf  JL)Hn  (V?  -L)exP(- 


where  m,  n  are  laser  modes,  H-^  and  Hn  are  Hermite  polynomials  of 
order  m  and  n. 


r  = 


(4) 


and 


and  rc  is  the  minimum  radius  attainable  (or  the  Gaussian  spot)  with  the 
coordinate  axis  defined  by 
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The  optical  arrangement  of  the  velocimeter  is  such  that  the  same 
lens  is  used  for  focusing  the  laser  beams  and  for  light  collection. 

Since  the  light  collecting  function  of  the  lens  is  characterized  by  a 
smaller  f/number  than  the  laser  focusing  function,  we  may  treat  the 
light  collecting  function  only  (since  its  depth  of  focus  is  considerably 
smaller  than  the  laser  beam  focal  depth).  Integrating  the  field  ampli¬ 
tude  (Ref.  2)  of  the  fundamental  laser  mode  yields  the  power  distribu¬ 
tion  as  a  function  of  radius  where  PQ  is  the  total  power  and  f(r)  is  some 
function  of  radius  r. 


P(r)  -  Po  (1  -  ef(r)) 


The  cross-sectional  area  of  the  beam  is 

A  =  nr  2 


and- with  the  use  of  Eq.  (4),  the  area  ratio  becomes 


Now, 


,  P(r) 
since  I  = 


where  I  is  the  intensity,  we  have 


laxial 

Imai 


(5) 


(6) 


(7) 


The  total  intensity  is  determined  by  letting  r-K°°  so  that  e^r^0  and 
Eq.  (7)  becomes 


^axi  al 

^max 


(8) 


This  equation  represents  the  axial  intensity  distribution  where 
Iaxial  is  the  intensity  at  an  axial  distance  of  z.  This  equation  was 
plotted  in  Fig.  4  for  the  argon  laser  of  A  =  4 8 80 A  for  a  variety  of 
f/numbers  where  rQ  was  calculated  with  the  conventional  diffraction 
limit  equation  for  minimum  focal  spot  diameter;  a  is  the  lens  diam¬ 
eter  and  L  its  focal  length. 
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Equation  (8)  can  also  be  simplified  by  rewriting  Eq.  (4)  in  the 
following  manner: 


which  is  identical  to  Eq.  (8)  so 

Ia*Ul  AoV 

CT  =  'r  '  (ID 


The  optical  configuration  of  the  LDV  is  determined  by  two  back-to- 
back  lenses  as  shown  in  the  sketch  below. 


If  a  particle  traverses  the  probe  volume  at  the  focal  point  of  lens  1,  the 
image  of  the  scattered  light  falls  on  the  aperture  preceding  the  phototube. 
If  the  aperture  is  made  equal  to  the  Gaussian  spot  size,  rQ,  then  the 
depth  of  the  probe  volume  is  given  by  Eq.  (8),  which  provides  the  small¬ 
est  possible  axial  dimension  that  can  be  achieved  for  a  particular 
f /number  lens.  If  the  aperture  size  is  larger  than  the  rQ,  a  greater  depth 
of  light  collection  will  be  effected.  Therefore,  it  will  be  necessary  to 
specify  the  axial  dimension  of  the  probe  volume  in  terms  of  the  aperture 
preceding  the  phototube. 

The  standard  thin  lens  equation  is  used  to  find  the  radius  of  a  spot 
on  the  aperture  plane  in  terms  of  light  scattering  particles  traversing  a 
distance  z  from  the  geometric  center  of  the  focal  point.  Let  the  focal 
plane  of  the  second  lens  be  the  position  where  the  aperture  will  later  be 
located.  Now,  starting  at  the  phototube  end  of  the  LDV,  we  can  deter¬ 
mine  the  shift  (z)  of  the  image  of  the  focal  spot  in  terms  of  the  z  distance 
from  the  focal  spot  in  the  object  plane. 
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Equation  (4)  governs  z'  in  terms  of  the  beam  radius  of  both  the 
Gaussian  spot  radius  and  any  arbitrary  radius  of  the  beam  at  some 
z  position.  The  parameters  r0  and  zQ  may  now  be  written  in  terms  of 
the  parameters  of  the  lens  configuration.  Thus 

c,  _  (£  +  |Sj)Li 

2  =  (t  +  |S{|)  -Lx  (12) 


Now,  Si' may  be  further  written  in  terms  of  the  lens  1  parameter 
and  the  z  distance  (that  distance  in  the  flow  field  that  a  particle  may 
traverse  off  the  geometric  center);  hence 


sr  = 


L](Li  +  z) 
z 


(13) 


When  Li  is  significantly  small  and  the  distance  is  large,  the  effec 
tive  beam  diameter  at  the  first  lens  must  be  taken  into  consideration 
such  that  ri  beam  =  divergence  times  total  image  distance;  so 


rl  beam 


riB  =  KSl  +  £)| 
Si 


(14) 


The  new  value  of  the  Gaussian  spot  in  terms  of  S2'is 


r„  = 


7i’rlB 


S2‘ 


(15) 


Substituting  into  Eq.  (4),  we  have 


since 


z'  =  |I<2  -  S2 


(16) 


(17) 


then 


(18) 


Equation  (18)  represents  the  radius  of  the  image  in  the  plane  of.the  aper¬ 
ture  (which  is  chosen  to  be  at  L2)  in  terms  of  the  two  lenses  of  the  LDV. 
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Prom  Eq.  (11),  the  axial  intensity  distribution  is  inversely  propor¬ 
tional  to  the  areas  of  two  separate  disuneters.  Hence 

/i»xi.i\  _  Aap\  2 

Vm  ax/  \r  )  (19) 

where  r^p  is  the  chosen  radius  of  the  aperture  preceding  the  phototube, 
and  r  the  radius  of  the  image  of  a  particle  traversing  the  probe  volume 
at  z.  Let  the  axial  dimension  of  the  probe  volume  be  defined  as  the  1/e^ 
of  Eq.  (20)  or  Eq.  (8);  then  z  may  be  found  from  Eq.  (17)  and  Eq.  (20), 
i.  e. , 


Again,  r^p  is  the  chosen  aperture  radius. 

The  total  dimension  is,  of  course,  2z,  since  z  was  defined  to  be 
zero  at  the  center  or  focal  spot  in  the  object  plane. 

Equation  (21)  may  be  approximated  by  assuming  a  linear  function 
between  the  new  focal  spot  point  and  the  aperture  plane,  i.  e. 

r  1 B  rAP 
Sz'  z' 

Writing  these  parameters  in  terms  of  the  input  parameters,  we 

have 


rAP  =  0.367  \llbl\ 
Si 


(22) 


where  rj  is  the  light  collecting  lens  radius,  Li  the  light  collecting  lens 
focal  point,  and  Si'  is  given  by  Eq.  (13).  Solving  Eq.  (22)  for  the  axial 
dimension  z  provides  a  good  first-order  approximation. 

i-i 

Z  =  _ 

[0.367  11-  -  lj  (23) 
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2.3.2  Theoretical  Analysis,  Radial  Dimension 

The  radial  dimension  of  the  data  volume  is  the  projected  aperture 
diameter  on  the  geometric  volume.  In  most  cases  the  lens  arrangement 
is  a  one-to-one  transfer  system  {no  magnification)  such  that  the  aperture 
diameter  is  the  diameter  of  the  ellipsoid. 

« 

It  was  previously  pointed  out,  and  will  be  shown,  that  an  even 
greater  reduction  in  the  z  dimension  could  be  effected  by  operating  the 
collecting  optics  off  axis  relative  to  the  axis  of  the  geometric  probe 
volume. 

The  geometric  probe  volume  and  the  data  probe  volume  are  expressed 
as  ellipsoids  at  the  1/e2  intensity  points.  For  clarity,  a  quarter  section 
of  each  probe  volume  is  shown  in  the  sketch  below  with  all  symbols  de¬ 
fined  as  in  the  above  equations.  * 


It  is  evident  from  the  sketch  that  the  axial  (z)  dimension  is  reduced 
by  2{z  cos  <£)  up  to  the  angle  <f>  =  a.  A  further  increase  in  <j>  (</>  >  a)  re¬ 
duces  the  axial  dimension  by 

P.V.L.  =  2(R  .cos  <&)  (24) 


1  / \x 

iOnly  the  cross  section  at  -y  =  0  will  be  considered  since  it  is 

assumed  that  the  volume  is  a  symmetric  ellipsoid  for  6  <  15  deg 
(see  Eq.  (2)). 
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It  is,  therefore,  only  necessary  to  determine  Rasa  function  of  the 
angle  <j> .  At  <f>  >  a,  R  varies  with  the  geometric  probe  volume.  The 
geometric  probe  volume  is  expressed  in  cylindrical  coordinates  as 


R  2cos  2<f>  +  R^sin^  _  ^ 

W  (-$  ■  (25> 


Solving  for  R,  we  have 


Substituting  the  terms  of  Az  and  Ay  from  Eq.  (2),  R  becomes 

R  =  b°  /2r 

yj  cos^<f)  sin2  0/2  +  cos 2  q/2 

With  the  use  of  the  "power  relation"  trigonometric  identity  and  multi¬ 
plying,  Eq.  (27)  reduces  to 


bo _ 

cos  6  cos  2<t>) 


(28) 


Equation  (28)  is  a  convenient  relationship  since  R  is  expressed  as  a 
function  of  the  laser  cross  beam  angle,  0  and  the  off-axis  angle,  $  . 


Combining  Eq.  (24)  with  Eq.  (28)  gives  the  probe  volume  length 
as 

P.V.L.  =  2b0  cos  * 

V  ‘/2(i  —  cos  6  cos  2 <j}) 

Let 

_ cos  ^ =  f(<p) 

yJW  -  cos  6  cos  2<£) 

then  the  length  of  the  data  probe  volume  becomes 

for  cj>>  a 
for  4>  <  a 


The  function,  f is  plotted  in  Fig.  5  versus  the  angle  <j)  for  a  variety 
of  laser  cross  beam  angles,  0.  To  use  the  curves  in  Fig.  5,  it  is  im¬ 
portant  to  calculate  a  to  determine  which  form  of  Eq.  (31)  to  use. 


P.V.L.  =  2b„f(<£) 
P.V.L.  =  2z  cos  0 


(29) 

(30) 

(31a) 

(32b) 
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At -the  Intersection  of  the  ellipsoid  and  circle,  described  by  the  ro¬ 
tation  of  the  data  probe  volume,  R  =  z,  and  $  =  a  {see  sketch).  In 
Eq.  (28),  setting  R  =  z  and  solving  for  cos  2a,  we  have 


cos  2 a 


zV2  -b* 
z2/2  cos  6 


(32) 


Rearranging  terms,  a  is  expressed  as 


a 


V2  cos-1 


p- 


-  2(b„/z)2 
cos  6 


(33) 


If  the  determination  of  a  requires  the  use  of  Eq.  (31b)  very  little 
reduction  is  realized  in  the  axial  dimension.  If  the  angle  a  is  very 
small,  small  off-axis  angles  will  produce  a  large  reduction  in  the  length 
of  the  probe  volume.  A  small  angle  <p  is  easily  produced  by  allowing  the 
two  laser  beams  to  fall  near  the  periphery  of  the  focusing/ collecting  lens 
(backs catter  units  only),  thereby  allowing  one  to  choose  a  very  small 
angle  <£ . 


When  a  is  desired  to  be  small,  or  negligible,  the  fraction  in  Eq.  (33) 
must  be  made  nearly  equal  to  unity,  such  that  cos  0  «  1  -  2(b0/z)2  or 

Y  =  >/'/,  (  1  -  cos  0)  <34) 

Hence,  one  can  choose  which  of  these  three  independent  parameters  to 
manipulate  to  establish  velocimeter  requirements. 

Equations  (31a)  and  (31b)  are  applicable  for  off-axis  angle  $,  up  to 
approximately  20  deg  since  the  value  R  is  not  the  outermost  extremity 
of  the  projected  data  probe  volume  or  the  geometric  probe  volume  axis. 
The  error  involved  for  angles  less  than  20  deg  is  assumed  to  be  negli¬ 
gible.  This  restriction  is  not  serious  since  the  trend  in  operational 
LDV  units  is  to  operate  in  the  backs  catter  mode  only.  In  this  case,  a 
20-deg  off-axis  operation  produces  such  low  scattering  intensities  that 
the  signal-to-noise  ratio  becomes  too  low  for  practical  velocimeter 
operation.  Another  reason  for  not  operating  at  angles  greater  than 
20-deg  is  demonstrated  in  Fig.  5.  For  angles  greater  than  15  deg,  the 
reduction  of  the  P.  V.  L.  becomes  extremely  small  regardless  of  the 
laser  cross  beam  angle. 
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SECTION  III 

EXPERIMENTAL  VERIFICATION 


3.1  BEAM  DIAMETER  DETERMINATION 

It  was  shown  (Fig.  1)  that  some  angular  relationship  governs  the 
reduction  of  the  axial  dimension  of  the  data  probe  volume.  The  experi¬ 
mental  curves  in  Fig.  1  could  not  be  predicted  theoretically.  An  error 
analysis  revealed  that  the  experimentally  determined  values  of  the  photo¬ 
tube  aperture  diameter,  2rA.Pi  and  the  laser  beam  diameter,  2b,  at  the 
focusing  lens  were  very  sensitive  to  the  axial  data  probe  volume  dimen¬ 
sions.  Heretofore,  2b  was  obtained  by  an  intuitive  guess,  since  accu¬ 
rate  measurements  of  2b  have  been  difficult,  if  not  impossible,  to  obtain. 
The  aperture  diameter,  2rAP»  was  assumed  to  be  circular,  which  gen¬ 
erally  was  not  the  case.  Also,  lint  or  other  particles  sometimes  blocked 
part  of  the  aperture  producing  smaller  effective  aperture  diameters. 
Therefore,  an  effort  was  directed  toward  the  accurate  measurement  of 
the  two  parameters.  The  aperture  diameter  measurement  problem  was 
solved  by  using  a  microscope  having  a  calibrated  gratical  in  the  eye 
piece.  This  permitted  direct  visual  verification  of  the  roundness  of  the 
aperture. 


The  first  attempt  to  determine  2b  was  to  count  the  number  of  cycles 
of  frequency  information  between  the  1/e2  points  in  the  center  of  the  geo¬ 
metric  probe  volume.  Since  N  (1/e2)  =  4/ 7r  D/2b,  where  D  is  the  two 
laser  beam  spacing  before  the  focusing  lens,  2b  may  be  calculated.  It 
was  found  that  this  technique  caused  a  10 -percent  error  in  the  deter¬ 
mination  of  2b,  hence  a  direct  measurement  was  undertaken. 


3.1.1  Wire  Technique 


A  method  used  for  the  direct  measurement  of  2b  was  to  pass  a  10— pi 
fiber  through  the  laser  beam  at  a  number  of  distances  from  the  laser, 
such  that  the  beam  divergence  could  be  determined.  A  fiber  was  then 
mounted  on  a  wheel  and  rotated  by  a  hysteresis -synchronous  motor. 
Figure  6a  shows  a  typical  laser  beam  scan  of  the  radial  diameter.  Par¬ 
ticular  notice  should  be  given  to  the  Gaussian  shape  of  this  fiber  scan. 
The  basic  assumption  of  the  fiber  passing  through  the  beam  is  that  the 
scattering  intensity  of  the  fiber  is  a  constant  multiple  times  the  Gaussian 
intensity  distribution  of  the  laser  beam.  To  prove  this,  we  must  show 
that  the  intensity  distribution  of  the  y-axis  (scan  on  y-axis)  is  Gaussian, 
i.  e. 


I(y)  =  fo  fj*  I/e  I(r)dx  dy 


(35) 
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when 

since  r^  =  x^  +  y^. 


or 


I(r)  =  ae-br2 

Substituting  I(r)  and  r^  into  I(y),  we  have 
I(y)  =  /  /  R  ae“b(x2  +  y2)dx  dv 

I(y)  =  ae-1>y2  *2  ae”bx  2  dx^  dy 


(36) 


The  integral  along  the  x-axis  would  be  a  constant  if  the  integration  were 
taken  from  zero  to  infinity.  But  the  integration  is  limited  to  the  l/e^ 
radius,  hence  Eq.  (36)  becomes 


l(y)  =  f(y)  (ae-by2) 


where  f(y)  is  the  weak  function  of  y.  Therefore,  (y)  is  not  a  Gaussian 
function  times  a  constant.  Since  f(y)  is  in  a  series  form  of 
(1  -  ay^  +  by^.  . . ),  the  error  is  small. 

A  series  of  data  points  taken  with  f(y)  assumed  to  be  constant  re¬ 
sulted  in  the  data  shown  in  Fig.  6b.  Several  photographs  were  taken 
(Fig.  6a)  at  several  distances  from  the  laser.  By  knowing  the  speed  of 
the  fiber  moving  through  the  laser  beam  and  the  sweep  speed  of  the 
oscilloscope,  the- 1/e  point  was  obtained  from  each  photograph.  (The  1/e 
point  was  chosen  since  the  background  radiation  affects  the  l/e^  intensity 
point  more  than  the  1/e  points.)  The  l/e2  point  was  calculated  by  assum¬ 
ing  a  Gaussian  intensity  distribution.  Figure  6b  shows  a  tremendous 
scatter  of  the  fiber  data  which  cannot  be  readily  explained  with  the 
assumption  of  f(y)  =  constant. 

3.1.2  Aperture  Technique  with  Phototube 

To  avoid  the  necessity  of  making  assumptions  as  in  section  3. 1.  1, 
an  alternate  technique  was  set  up  to  measure  the  laser  beam  diameter 
versus  distance.  A  10-A*  aperture  preceding  a  phototube  was  mounted 
on  a  microdensitometer  traversing  table.  In  addition  to  the  horizontal 
traverse,  vertical  adjustment  was  also  incorporated  such  that  each 
position  in  the  beam  could  be  optimized  to  ensure  the  scan  was  on  the 
diameter  of  the  beam. 

The  plots  of  the  data  indicated  a  s  atur  at  ion  /  f  atigue  effect  after  scan¬ 
ning  through  the  maximum  intensity.  Figure  4a  indicates  this  effect  very 
clearly.  To  retain  these  data,  a  computer  curve -fitting  technique 
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(Appendix  EE)  was  applied  only  to  the  data  on  the  fatigued  side  of  the 
curve  and  to  points  approximately  10  percent  over  the  maximum.  The 
resulting  l/e^  points  are  plotted  in  Fig.  6b. 

3.1.3  Aperture  Technique  with  Diode 

It  was  believed  that  the  light  intensity  passing  through  the  10-M 
aperture  would  be  too  weak  for  conventional  amplitude  measuring  de¬ 
vices,  hence  a  phototube  was  incorporated  in  the  measurement.  Be¬ 
cause  of  the  fatigue /saturation  problem,  one  datum  point  was  taken  with 
a  photodiode.  This  technique  requires  a  more  elaborate  experimental 
setup  since  two  separate  traverse  units  are  required.  The  procedure 
followed  was  identical  to  that  of  the  aperture  phototube  technique  de¬ 
scribed  above.  Figure'  6c  shows  the  data  acquired  7  ft  from  the  laser. 
Included  in  this  curve  is  the  computer -generated  curve  fit  which  clearly 
shows  the  superiority  of  the  technique. 

Based  on  the  diode  data  only,  the  following  empirical  relationship 
giving  the  beam  diameter  versus  distance  was  derived  for  the  Spectra 
Physics  124  laser 

(2b±0.005  =  0.11  x  +  1-Ds.P.  124  (37) 

where  x  is  the  distance  from  the  laser  in  feet,  and  2b  is  the  beam  diam¬ 
eter  in  millimeters. 

The  argon  laser.  Spectra  Physics  140,  has  been  the  primary  laser 
used  for  dual  scatter  LDV  work.  A  similar  equation  for  beam  diameter 
versus  distance  would  also  be  required.  During  the  course  of  the  experi¬ 
ment,  considerable  difficulty  was  encountered  in  adjusting  the  laser  for 
TEM00  mode  operation.  Therefore,  an  automatic  diode  scan  technique 
was  established  for  rapid  determination  of  the  laser  mode  structure.  To 
ensure  sufficient  accuracy  in  this  determination,  several  scans  were 
made  at  19.  5  ft.  A  typical  scan  is  shown  in  Fig.  7.  The  resulting  em¬ 
pirical  relationship  was  determined  to  be 

(2b=0.05  =  0.177  x  +  1.6)s.p.  140  (38) 

where  x  is  in  feet  and  2b  is  in  millimeters. 


3.2  DATA  PROBE  VOLUME  DETERMINATION 

Verification  of  Eq.  (31)  was  effected  primarily  with  a  backscatter 
LDV  unit.  A  schematic  of  a  typical  off-axis  lens  arrangement  is  shown 
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in  Fig.  8.  The  on-axis  experiments  also  used  backscatter  units  since 
this  eliminated  the  necessity  of  blocking  the  two  cross  beams  so  they 
would  not  enter  the  phototube.  The  off-axis  experiments  were  all  con¬ 
ducted  with  apertures  slightly  less  than  the  Ay  dimension  of  the  geomet¬ 
ric  probe  volume.  This  allowed  direct  comparisons  between  off-axis, 
large  aperture  operation  and  on- axis,  small  aperture  operation.  It  is 
desirable  to  operate  the  velocimeter  with  large  aperture  because  of  in¬ 
creased  signal  amplitude,  increased  signal -to-noise  ratio,  more  cycles 
in  the  information  burst,  and  greater  vibration  tolerances.  A  tabulation 
of  the  data  of  the  on- axis  and  off-axis  experiments  is  included  in  Table  I. 
Again,  0  is  the  off-axis  angle.  The  length  of  the  theoretical  probe  vol¬ 
ume  was  calculated  from  the  parameters  of  each  experiment,  and  the 
experimental  probe  volume  length  was  measured  directly. as  described 
in  Ref.  1.  Excellent  agreement  can  be  seen  in  P.  V.  L.  experiment  and 
P.  V.  L.  theory  ^  was  ^at  ma3or  error  of  these  experiments 
was  manifested  in  the  angle  0  because  it  was  extremely  difficult  to  de¬ 
termine. 


SECTION  IV 

SUMMARY  AND  CONCLUSIONS 


It  was  shown  in  Fig.  1  that  the  axial  dimension  of  the  probe  volume 
could  be  considerably  reduced  by  operating  the  velocimeter  off  axis.  All 
of  the  experiments  were  conducted  with  approximately  the  same  geo¬ 
metric  probe  volume  length  of  0.  8  in.  This  allows  one  to  compare  the 
results  in  Fig.  1  in  relation  to  the  relatively  long  probe  volume  length.  A 
vital  point  must  be  remembered  in  the  design  of  a  velocimeter:  velocity 
information  is  present  in  the  entire  geometric  probe  volume  length,  but 
with  a  lens/ aperture  combination  or  off-axis  operation  only  a  portion  of 
the  total  available  information  is  used.  In  many  cases,  operating  off 
axis  with  reduced  scattering  intensity  produces  greater  signal  intensity 
and  more  cycles  of  the  desired  information  than  operating  on  axis  with 
very  small  apertures.  In  general,  smaller  probe  volume  lengths  are 
obtainable  (a  very  desirable  feature  for  velocity  resolution)  with  small 
off-axis  operation  as  shown  in  Fig.  1  and  elucidated  further  in  Fig.  5. 

Figure  5  provides  a  good  illustration  of  the  "best  possible"  probe 
volume  reduction  for  any  off-axis  angle  and  any  beam  spacing,  0, 
desired  for  velocimeter  design.  It  is  shown  that  a  beam  spacing  of 
greater  than  15  deg  would  produce  practically  no  reduction  in  the  length 
of  the  probe  volume  by  operating  off  axis.  On  the  other  hand,  a  tremen¬ 
dous  reduction  is  seen  for  a  beam  spacing  of,  say,  0.  5  deg  and  very 
small  off-axis  angles. 
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TABLE  I 

RESULTS  OF  PROBE  VOLUME  EXPERIMENTS 


n 

2  rAP 

9,  deg 

2  b,  nun 

m 

*(<*>) 

P.V.L.  Theoretical 

P.V.L.  Experimental 

0  . 

0.01375 

in. 

— 

— 

(0.377)  in. 

Eq.  (23) 

0.36  in. 

0 

165  p. 

— 

— 

(0.178)  in. 

Eq.  (23) 

0.188  in. 

8 

0.018 

in. 

2.67 

1.6 

0.0151 

in. 

7.70 

(0.107)  in. 
Aperture  too  large 

0.247  in. 

8.5 

465  ii 

2.633 

1.43 

0.426 

mm 

6.61 

(2.817  nun) 

0.011  in. 

(2.82  aim) 

0.01  in. 

10 

0,018 

in. 

2.62 

1.50 

0.02346 

in. 

5.50 

0.129  in. 

0.127  in. 

13 

0.018 

in. 

2.70 

1.3 

0.0186  ± 
0.002  in. 

3.75  ±  0.3 

0.0697  ± 

0.0131  in. 

0.073  ±  0.004 
in. 

14 

465  (i 

2.50 

1.14 

0.02120  ± 
0.001  in. 

3.839 

0.0814  + 

0.008  in. 

0.0919  ±  0.016 
in. 
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An  interesting  fact  is  borne  out  in  Eq.  (31),  which  shows  that  the 
probe  volume  length  is  independent  of  the  parameters  of  the  collector 
optics.  Since  a  is  affected  directly,  care  must  be  taken  in  velocimeter 
design  not  to  exceed  the  range  of  <£.  Conversely,  for  angles  less  than  a, 
P.  V.  L.  is  independent  of  the  laser  cross  beam  parameters. 

The  resulting  probe  volume  in  off-axis  operation  is  obviously  not 
symmetrical  compared  to  the  on- axis  operation,  but  it  is  felt  that  this 
is  of  little  or  no  consequence  since  the  important  factor  in  reducing 
the  probe  volume  length  is  velocity  resolution,  not  symmetry. 

The  probe  volume  length  is  greatly  affected  by  the  electronic  read¬ 
out  unit  such  that  the  P.  V.  L.  may  be  further  reduced  or  in  some  cases 
even  increased.  The  only  effect  the  electronics  have  on  the  probe  vol¬ 
ume  is  in  the  form  of  discrimination  of,  for  example,  amplitude  ratios 
of  the  information  pulse  and  frequency  deviation.  The  most  effective 
method  is  to  use  a  variable  electronic  threshold  to  reject  all  signal 
amplitudes  below  a  preset  level  for  the  axial  intensity  distribution  of  the 
probe  volume.  As  much  as  90-percent  reduction  could  be  effected  by 
this  method. 

The  acceleration  of  the  light  scattering  particles  passing  through 
the  data  probe  volume  must  be  negligible  for  this  method  to  be  effective. 
Other  effects  of  electronic  probe  volume  reduction  are  generally  unique 
to  the  type  of  readout  used  in  the  measurement.  Two  specific  examples 
of  readout  effects  on  probe  volume  are  worth  mentioning:  the  spectrum 
analyzer  and  the  zero  crossing. 

The  spectrum  analyzer  is  a  common  technique  for  providing  the 
LDV  velocity  information.  These  analyzers  display  the  power  spectrum 
of  the  frequency  (velocity)  content  of  the  data  probe  volume.  Assuming 
negligible  acceleration  in  the  flow,  seeded  with  a  size  distribution  having 
a  narrow  variance  (not  too  difficult  to  obtain  in  practice),  discrimination 
is  effected  by  recording  the  maximum  intensity  of  the  distribution  func¬ 
tion,  since  particles  moving  at  different  velocities  relative  to  geometric 
center  of  the  data  probe  volume  produce  lower  amplitude  signals.  This 
type  of  amplitude  discrimination  may  produce  extremely  small  probe 
volume  lengths,  such  as  experiments  concerned  with  measuring  veloci¬ 
ties  through  shock  waves  where  the  sudden  "break"  in  the  velocity  pro¬ 
file  could  easily  be  determined.  It  was  experimentally  determined  that 
a  probe  volume  movement  of  0.001  in.  produced  this  "break"  in  the 
velocity  profile.  Hence,  for  comparison. 

Geometric  probe  volume  length  0.  9  in. 

Data  probe  volume  length  0.  2  in. 

Effective  probe  volume  length  0.  001  in. 
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Another  experiment  using  a  zero  crossing  readout  was  used  with 
the  same  flow  and  LDV  unit.  The  zero  crossing  unit  reads  the  time  it 
takes  a  frequency  burst  to  cross  zero  level,  say,  ten  times.  Obviously, 
this  type  of  readout  (in  its  basic  form)  cannot  be  very  sensitive  to  ampli¬ 
tude  discrimination;  regardless  of  the  amplitude  of  the  information  sig¬ 
nal,  a  "perfect"  zero  crossing  readout  will  provide  the  velocity  informa¬ 
tion.  This  means  that  regardless  of  the  size  of  the  data  probe  volume, 
the  readout  will  "read"  all  velocities  in  the  geometric  probe  volume. 

This  was  shown  experimentally  to  be  a  case  with  the  same  geometric 
probe  volume  and  data  probe,  volume  as  in  the  above  example 
(G.  P.  V.  L.  =0.9  in;  D.  P.  V.  L.  =  0.  2  in).  The  flow  was  a  0.  5-in.-diam 
free  jet  having  a  nearly  constant  velocity  in  the  core  of  the  flow.  The 
zero  crossing  readout  produced  the  velocity  distribution  shown  in  Fig.  9 
with  the  corresponding  spectrum  analyzer  trace.  The  figure  shows  the 
mean  velocity  in  the  center  of  the  data  probe  volume  along  with  a  portion 
of  velocities  in  the  geometric  probe  volume.  Since  the  geometric  probe 
volume  is  much  greater  than  the  flow  field,  all  velocities,  from  zero  to 
maximum  centerline  velocity,  should  be  present.  The  decreasing  distri¬ 
bution  may  be  explained  by  the  fact  that  complete  independence  from 
amplitude  discrimination  is  not  feasible  in  "real"  electronic  readout 
units.  Hence,  the  shape  of  the  velocity  distribution  may  be  explained  by 
the  near-perfect  readout  analogy.  Needless  to  say,  a  threshhold  built 
into  this  readout  would  again  recover  the  amplitude  discrimination  for 
obtaining  very  small  effective  probe  volumes. 
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APPENDIXES 


I.  ILLUSTRATIONS 
II.  GAUSSIAN  CURVE  FITTING 


Normalized  Amplitude 


Fig.  1  Experimental  Data  of  On-Axis  and  Off-Axis  P.  V.  L. 
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Fig.  2  Defining  Parameters  of  Beam  Crossover  Region 
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light  shield. 


Curtain 
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View  A-A 
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by  geometry  of  the  lens 
configuration  as  shown. 


Fig.  3  Light  Shield  for  Reducing  Background  Light  Effects 
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Probe  Volume  Length  (PVL)  is  given  by 
PVL  -  2b0  f(0>,  (0  >  a) 

and 

PVL  -  2z  cos  0,  (0  <  o) 


The  percentage  of  Probe  Volume  Reduction  (A)  for 
specified  0  is  given  by 


A%  ■  100 


MtH.> 


and 


A%  ■  100  [l  -  cos  0J0  <  a 
where 

0  ■  Off-axis  angle,  deg;  0  >  0 


o-  1  cos"1 


l-2{2 
z 


cos  e 


f<0> 


cos  0 


VF 


cos  8  cos  20) 


z « A  Ion-axis  data  probe  volume) 


5  10 

Off-Axis  Angle,  0,  deg 

Fig.  5  Plot  of  the  Function  f (0)  versus  4> 
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Typical  Photograph  of  Fiber  Passing 
through  Laser  Beam 


Typical  Data  Curve  from 
Aperture  with  Phototube 
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Radial  Position 


a.  Representative  Data  of  Laser  Beam  Intensity  versus  Radial  Position 

Fig.  6  Typical  Data  Plots 
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b.  .Plot  of  Beam  Diameter  versus  Distance  for  SP-124 
Fig.  6  Continued 
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Radial  Position 


c.  Intensity  Plot  of  Diode  Data  with  Computer  Curve  Fit 
Fig.  6  Concluded 


(Scan  Speed)(l/2  Width)  x  1.428 
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Fig.  7  Automatic  Diode  Scan  of  Laser  Beam  Diameter  for  SP-140 
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a  +  b  =  x  for  Eq.  (37) 

2  rAp  =  Aperture  Diameter 


Fig.  8  Typical  Setup  for  Off-Axis  Experiments 
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Number  of  Occurrences 


Corresponding  Spectrum 
Analyzer  Trace 


Time 


4.77  4.35  4.00  3.71  3.43  3.23  3.03  2.86  2.70  2.56 


5.00  4.54  4.17  3.85  3.57  3.33  3.13  2.94  2.78  2.63  2.50 

Frequency 


Fig.  9  Zero  Crossing  Electronic  Readout 


AEDC-TR-71-85 


AEDC-TR-71  -85 


APPENDIX  II 

GAUSSIAN  CURVE  FITTING 


Assume  that  experiment  provides  a  table  of  the  quantities  I  and  x, 
and  that  I  is  known  to  be  a  Gaussian  function  .of  x. 


,  ,  — 2(x  -  x0)2 

I  =  I0e  - 

b2 


(II- 1) 


where  IQ  is  the  maximum  value  of  I,  xQ  is  the  origin  about  which  the 
x  data  are  taken  (only  in  rare  cases  would  x0  =  0),  and  b  is  a  constant. 
The  constant  b  is  useful  because  I(x  =  b)  =  1 0/e2  and  therefore  2b  is  a 
convenient  measure  of  the  "diameter"  of  the  Gaussian  function. 


This  appendix  describes  how  one  can  determine  the  value  of  b  from 
a  table  of  I  versus  x.  The  method  of  least  squares  is  employed. 


Most  computing  facilities  have  a  standard  program  for  performing 
least-squares  curve  fits  to  function  of  the  form 

y  =  a0  ax  x  &2  x2  +  •••  +  an  xn  ( II  —  2) 


in  which 


and 


y  =  In  I 


aQ  —  In  I0 


al 


4l« 

b2 


2X02 

b2 


(n-3) 

( II  -  4 ) 
( II  -  5 ) 

(II  -  6 ) 


A  table  of  I  versus  x  containing  N  entries  can  be  written  in  the  matrix 
form 


or  as 


yi  " 

1  Xl  XI 2 

y2 

1  x2  x22 

"ao" 

• 

- 

■  •  ■ 

al 

•  •  • 

•  •  • 

a2 

_>'n_ 

1  xn  xn 

(H-7) 


Y  =  XA 


(II -8) 
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The  matrix  A  is  found  to  be 


A  =  (X  X)“l  X  Y 

(II  -  9 ) 

The  variance  in  y  is  given  by 

„  2  _  YY  -  AXY 
y  N-l 

(11-10) 

and  the  variance-covariance  matrix  (the  diagonal  elements 
variances  in  the  constants  aQ,  aj,  and  a2>  is  given  by 

are  the 

V  =  oy2  (X  X)“l 

(II- 11) 

In  order  to  apply  the  above  matrix  equations,  one  must  prepare  a 
table  of  In  I  versus  x  from  the  original  table  of  I  versus  x,  and  then 
write  the  new  table  in  matrix  notation.  The  matrix  manipulations  are 
then  done;  one  finds  the  Gaussian  diameter  from  Eqs.  ( II - 9 )  and  ( II - 6 ) . 

2b  =/=J 

V  o2 

(11-12) 

The  standard  deviation  in  the  diameter  is  given  by 

a  (2b)  a  (a2) 

a3/2 

-b3  (7(02) 

0 

(11-13) 

where  0^2)^  is  the  appropriate  element  of  the  variance-covariance 
matrix  V. 


33 


UNCLASSIFIED 

Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 

f Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  la  dossMed} 


3  REPORT  TITLE 

LASER  DOPPLER  VELOCIMETER  DUAL  SCATTER  PROBE  VOLUME 


4  descriptive  NOTES  (Type  of  report  end  Inclusive  datea) 

Final  Report  -  July  through  October  1970 

5  autmoRiSj  (Flral  name,  middle  initial,  feat  name) 

W.  H.  Goethert,  ARO,  Inc. 


10  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13  ABS’RAC — 

A  further  clarification  of  the  volume  from  which  data  originate 
in  the  dual  scatter  laser  Doppler  velocimeter  is  presented.  There 
are  several  factors  that  specify  the  primary  dimensions  of  the  probe 
volume  which  in  some  cases  are  totally  specified  by  the  collector 
optics.  It  was  found  experimentally  that  the  lens  aperture  of  the 
light  collecting  optics  has  a  primary  effect  on  the  volume  from 
which  data  originate.  An  equation  is  derived  for  reducing  this 
data  volume  for  both  the  on-axis  and  off-axis  cases.  Experimental 
verification  is  presented  and  compared  to  the  derived  equation. 


DD  ,'“'".,1473 


11  SUPPLEMENTARY  NOTES 

Available  in  DDC 


12  SPONSORING  MILITARY  AC  Tl  VITY 

Arnold  Engineering  Development 
Center,  Air  Force  Systems  Command, 
Arnold  AF  Station,  Tenn.  37389 


ft  REPOR  T  O  A  TE 

July  1971 

Bm.  CONTRACT  OR  GRANT  NO 

F40600-72-C-0003 

b.  PROJEC  T  NO  4344 

Task  32 

d  Program  Element  65701F 


7a.  TOTAL  NO  OP  PAGES  7b.  NO  OF  REFS 

39  4 

9*.  ORIGINATOR'S  REPORT  NUMBER(S) 

AEDC-TR-71-85 


9b.  OTHER  REPORT  NO(5>  ( Any  other  numbara  that  may  bm  aaatgned 
•  this  report) 

AR0-0MD-TR-71-24 


I  ORIGINATING  ACTIVITY  (COfpOtMI.  author) 

Arnold  Engineering  Development  Center 

ARO,  Inc.,  Operating  Contractor 

Arnold  Air  Force  Station,  Tennessee  37389 


2a.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

2b.  GROUP 

N/A 


UNCLASSIFIED 

Security  Classification 


UNCLASSIFIED 


